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We present a Green's function based framework for modeling the scanning tunneling spectrum 
from the normal as well as the superconducting state of complex materials where the nature of 
the tunneling process— i.e. the effect of the tunneling 'matrix element', is properly taken into 
account. The formalism is applied to the case of optimally doped Bi2Sr2CaCu2 0s+5 (Bi2212) high- 
Tc superconductor using a large tight-binding basis set of electron and hole orbitals. The results 
show clearly that the spectrum is modified strongly by the effects of the tunneling matrix element 
and that it is not a simple replica of the local density of states (LDOS) of the Cu-d x z_ y 2 orbitals 
with other orbitals playing a key role in shaping the spectra. We show how the spectrum can be 
decomposed usefully in terms of tunneling 'channels' or paths through which the current flows from 
various orbitals in the system to the scanning tip. Such an analysis reveals symmetry forbidden and 
symmetry enhanced paths between the tip and the cuprate layers. Significant contributions arise 
from not only the CuCb layer closest to the tip, but also from the second CuC>2 layer. The spectrum 
also contains a longer range background reflecting the non-local nature of the underlying Bloch 
states. In the superconducting state, coherence peaks are found to be dominated by the anomalous 
components of Green's function. 



PACS numbers: 68.37.Ef 71.20.-b 74.50. +r 74.72.-h 

I. INTRODUCTION 

High resolution scanning tunneling spectroscopy (STS) 
together with other highly resolved spectroscopies such 
as angle resolved photoemission (ARPES), is making it 
possible to obtain a comprehensive mapping of the elec- 
tronic spectrum of the high-temperature superconduc- 
tors (HTSs) in both real and reciprocal space over a 
wide range of dopings and temperatures. These stud- 
ies are providing insight into the rich phase diagrams of 
the HTSs, and are leading thus to an understanding of 
the 'missing links' for developing a definitive theory of 
how high superconducting transition temperatures arise 
in these unconventional materials. In STS experiments, 
the focus to date has been on hole doped cuprates, es- 
pecially on Bi2Sr2CaCu208+5 (Bi2212), which has been 
the subject of an overwhelming amount of experimental 
work, see, e.g, Refs. HH6J Bi2212 is a typical cuprate 
material, which is an antiferromagetic insulator in the 
strongly underdoped regime, but exhibits a supercon- 
ducting phase over a wide range of hole doping. 



STS can be applied to a substantial part of the doping 
and temperature spanned phase space of HTS materials. 
The superconducting (SC) phase is observed around op- 
timal hole doping (OP), while the pseudogap (PG) phase 
is found within the underdoped regime (UD). As a prac- 
tical limitation, STS requires a conducting sample, but 
the deeply underdoped regime is insulating and hence 
unreachable by STS. However, under experimental condi- 
tions the samples are not homogeneously doped. Rather, 
there is a strong spatial variation in doping, which makes 
observation of a continuum from the PG to the SC phase 
possible within one sample. Although these spatial vari- 
ations in STS generally appear irregular, quite recently a 
more ordered coexistence of PG and SC phases has been 
observed^. 

The physics of the cuprates is dominated by the 
cuprate layers, which are usually not exposed to the tip 
of the apparatus. For example, in Bi2212, the quasipar- 
ticle tunneling takes place through insulating BiO and 
SrO layers. The conventional interpretation of the spec- 
tra is based on the assumption that the STS spectrum 
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is directly proportional to the LDOS of the C11O2 layer, 
especially the LDOS of the d x i_ y i orbitals, thus neglect- 
ing the effects of the tunneling process in modifying the 
spectrum in the presence of the insulating overlayers and 
multiple orbitals. The motivation for this simplification 
is an attempt to reduce the quasiparticle structure to few 
band models, which are amenable to theoretical treat- 
ment of strong correlation effects in the presence of super- 
conducting and antiferromagnetic order. Notably, there 
have been attempts to take the effect of the overlayers 
into account by assuming a 'tunneling matrix element' 
or a 'filter function' 1 ^ 9 !. 

With this background, our recent work on STS 10 
of Bi2212 provides a significant advance in realistic 
material-specific modeling of the STS spectrum. We in- 
voke a Green's function approach where a large number 
of orbitals is included, and all tunneling paths to the 
tip in the semi-infinite solid are taken into account. We 
showed clearly that instead of being a simple reflection 
of LDOS of the Qw-d x i_ y i orbitals, the STS signal repre- 
sents a very complex mapping of the electronic structure 
of the system. 

In this study we extend our approach by decomposing 
the tunneling current in terms of regular and anomalous 
matrix elements of the spectral function in an atomic or- 
bital basis. As in Ref. [TDJ we concentrate on Bi2212 
as the canonical HTS material. We start by reformu- 
lating the well-established methods to model tunneling 
current in nanostructures into a more transparent form 
for interpreting tunneling in the superconducting state. 
Our derivation is based on the conventional Todorov- 
p en( ^ r ^iili2j approach (TP), which is closely related to 
the more common Tersoff-HamanrPI method (TH). TP 
and TH methods both employ a calculation of the LDOS, 
but TP is more naturally written in terms of Green's 
functions. We will show, in fact, that TP decomposes 
into matrix elements of the spectral function, giving very 
detailed information concerning the origin of various fea- 
tures in the tunneling spectrum. We thus demonstrate 
how the contribution of different atomic orbitals to the 
total current can be extracted from the calculations. Our 
spectral decomposition also naturally distinguishes be- 
tween the electron and hole nature of the quasiparticles 
in the superconducting state. In addition, it leads to 
a multiband generalization of filtering function by Mar- 
tin et aP and a clarification of selection rules governing 
tunneling through filtering layers. This information is 
important, e.g., in determining how a dopant or impu- 
rity atom alters the spectrum, and how the effect of such 
a perturbation is seen in real space. 

In order to gain a handle on the effects of filtering 
layers, we derive a consistent form of a filter function 
through Green's function manipulations. This rigorous 
form for the filtering effects is useful for determining the 
relation between the tunnel current and the LDOS of the 
Cu02 layers. We show that this relation is nontrivial in 
that some channels are 'first-order forbidden'. Thus our 
new approach shows that no direct regular signal from 



d x 2_ y 2 orbitals of the Cu directly below the STM tip 
reaches the microscope. Instead the d x i_ y i orbitals of 
the four neighboring Cu atoms give a major contribu- 
tion to the tunneling signal. Although we concentrate 
on pristine systems in the present work, the results have 
important implications for inhomogeneous situations - 
e.g., the relationship between the observed features in 
the spectrum of an impurity atom and the underlying 
LDOS. This decomposition also allows treatment of the 
regular and anomalous propagation of quasiparticles in 
a superconductor, and on this basis we show that the 
coherence peaks result from the anomalous electron-hole 
propagation. 

The paper is organized as follows. The model for the 
geometrical structure and the electronic structure is in- 
troduced in Sections II. A and B, respectively. The meth- 
ods to calculate the Green's function in the normal and 
the superconducting state are derived in Sections II. B 
and C, respectively. The Todorov-Pendry equation for 
the tunneling current is decomposed into regular and 
anomalous terms to show not only the proper form of the 
matrix element but also the partial current terms for any 
chosen orbital in Section II. D. The formalism is applied 
to discuss STM topographic maps in Section III. A, and 
the STS spectrum of Bi2212 in Section III.B. The spec- 
trum is then analyzed in terms of tunneling matrix ele- 
ments and partial currents in Section III.C. Further com- 
ments on symmetry analysis are made in Section IV. A., 
and remarks on electron extraction/injection are made in 
Section IV. B. Finally, conclusions are drawn and future 
applications sketched in Section V. Relevant technical 
details of the form of boson-electron coupling assumed 
in the tunnel spectra and of the superconducting state 
calculations are given in the two appendices. 

II. DESCRIPTION OF THE MODEL 

Our theoretical framework involves three distinct 
steps. First, we choose a three-dimensional geometri- 
cal model of atoms with a sufficiently large simulation 
cell with periodic boundary conditions in the horizontal 
directions to treat a semi-infinite solid surface. Second, 
we attach a basis set of atomic orbitals to each atom. At 
this stage, the one-particle Hamiltonian is constructed 
and the corresponding Green's function tensor is formed. 
Third, we apply our Green's function formalism to eval- 
uate the tunneling current. The technical details of these 
three steps are outlined in the following three subsections. 

A. Sample geometry 

We model the Bi2212 sample as a slab of seven layers 14 
in which the topmost layer is BiO, followed by layers 
of SrO, Cu02, Ca, Cu02, SrO, and BiO, as shown in 
Fig. [TJa). The tunneling computations are based on a 
2V2 x 2V2 real space supercell consisting of 8 primitive 
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FIG. 1: (color online) (a) Side view of the simulation cell 
used to compute the tunneling spectrum of Bi2212. Tunnel- 
ing signal from the conducting CuCb layers reaches the tip 
after passing through the filtering layers of SrO and BiO. (b) 
Cuprate layer showing the supercell consisting of eight prim- 
itive cells, (c) Top view of the surface showing the arrange- 
ment of various atoms. 



surface cells with a total of 120 atoms (see Fig. [TJb)) . 
The coordinates are taken from the tetragonal crystal 
structure of Ref. 15. For STS simulations, the STM 
tip is modeled as an orbital with an s-wave symmetry at 
the assumed position of the apex of the tip. This tip is 
allowed to scan across the substrate for generating the 
topographic maps such as those in Fig. [5j or held fixed 
on top of a surface Bi atom for the computed spectra 
presented for example in Fig. [6] 



B. Construction of the uncorrelated normal state 
Hamiltonian 

In order to construct a realistic framework capable of 
describing the tunneling spectrum of the normal as well 
as the superconducting state of the cuprates, we start 
with the normal state Hamiltonian for the semi-infinite 
solid in the form 



H, = 



a(3cr 



[e a ci a c aa + V a/3 cl a cp a ] 



(1) 



which describes a system of tight-binding orbitals created 
(or annihilated) via the real-space operators (or c aa ). 
Here a is a composite index denoting both the type of 
orbital (e.g. Cu-d x 2_ y 2) and the site on which this or- 
bital is placed, and a is the spin index. e a is the on-site 
energy of the a th orbital, a and (3 orbitals interact with 
each other through the potential V a p to create the energy 
eigenstates of the entire system. 

The specific electron and hole orbital sets used for var- 
ious atoms are: (s,p x ,p y ,p z ) for Bi, Ca and O; s for Sr; 



and (4s, d 3z 2_ r 2,d xy , d xz , d yz , d x 2_ y 2) for Cu atoms. This 
yields 58 electron or hole orbitals in a primitive cell and 
a total of 2 x 464 orbitals in the 2\/2 x 2^/2 simulation 
supercell. The number of k-points used in the compu- 
tations depends on whether we do band calculations or 
solve the Green's function. For band calculations, we use 
a dense set of k- values to produce smooth bands for di- 
rections T^M^X^Tas seen for example in Fig. 
|2j In the case of Green's function calculations, we use 
Nk = 256 k-points for the supercell Brillouin zone. This 
corresponds to 8 x 256 = 2048 k-points for a primitive 
cell. 
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TABLE I: Slater-Koster prefactors, v a p m , and onsite energies 
s a . The Vcxp m are used to construct the Hamiltonian overlap 
matrix elements V a p as described in Ref.™ 

The Slater-Koster formalismP^L^ is used to fix the an- 
gular dependence of the tight binding overlap integrals. 
The onsite energies and the prefactors are fitted to the 
LDA band structure of Bi2212 that underlies for example 
the extensive angle-resolved photointensity computations 
of Refs. 19-24. In Table I, we show the specific values of 
the v a /3 m prefactors used for computing the Slater-Koster 
hopping integrals. Notably, we have shifted the bottom 
of the BiO conduction band to agree with experiments, 
which do not observe the Bi-bands at least within leV 
above the Fermi-level. This choice is also supported by 
calculations of Ref. [25] which show the sensitivity of the 
position of the Bi-band with respect to impurities and 
doping. The absence of the bottom of the BiO band in 
the STS spectra may also be due to a voltage gradient 
across the insulating filter layers (BiO and SrO layers) 
when applying a bias voltage between the tip and the 
sample. If so, the absolute value of the voltage within 
these layers is less than the bias voltage T4> and thus the 
apparatus would need to apply a bias which would be 
significantly larger than Vb to locally see states that are 
strictly at Ep + eV^. 

The tight-binding parameters of the normal state 
Hamiltonian of Eq. ([!]) produce the detailed band struc- 
ture of Bi2212 shown in Fig. [2] While the tight-binding 
band structure is in reasonable agreement with the LDA 
band structure of Ref. [25j in order to carry out spec- 
troscopic computations, one must additionally make sure 
that the underlying wavefunctions are described correctly 
including their symmetries. Our procedure based on the 
the use of Koster- Slater matrix elements not only fits the 
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band stuctures, but the symmetries and phases of the 
associated wavefunctions are also described correctly. 
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FIG. 2: (color online) (a)-(d): Normal state band structure 
of Bi2212 for the tight-binding Hamiltonian and from first- 
principles LDA computations. Weights of Cu-d x 2_ y 2 and Cu- 
d z 2 contribution to the bands are shown using a color scale 
where red denotes high and blue low values (see colorbar). 
Note that the tight-binding calculations are done for a slab, so 
that the the tight-binding bands do not display the splitting 
of Bi-0 bands seen in LDA results. The LDA bands have 
been calculated using Virtual Crystal Approximation (VCA) 
with 24% Pb doping to set the bottom of the BiO band, (e) 
and (f): Quasiparticle band structure in the superconducting 
state based on the Hamiltonian of Eq. (|3| is shown in (e). 
Panel (f) zooms in on the gap region of (e) which is shaded 
grey. Electron character of quasiparticles is shown in red and 
the hole character in blue. Notice, that the quasiparticles 
differ significantly from being electrons or holes only in the 
close neighborhood of the superconducting gap around the 
M-point. 

Figs. [2] (a)-(b) show the normal state tight-binding 
band structure based on our 58 orbital Hamiltonian of 
Eq. ([T]). The main cuprate bands, with predominantly 
Cu-d x 2_ y 2 character, are seen in panels (a) and (b) to fol- 
low the corresponding LDA calculations in panels (c) and 
(d). Note that in our tight-binding modeling, we have 



adjusted the positions and bilayer splitting of the two 
van Hove singularities (VHSs) to approximately match 
the experimental photoemission and STS findings for 
the optimal doping (OP) region with hole concentration 
p « 0.16 (Ref. I26|27p . In addition to Cu-d x 2_ y 2, Cu-d z 2 
is seen in panels (b) and (d) to give a significant spec- 
tral weight to this band, especially at energies below the 
Fermi- level. The complicated 'spaghetti' region has large 
contributions from the d z 2 of Cu and horizontal p x (p y ) 
orbitals of the oxygens within the cuprate layer as well as 
the vertical p z orbital of the apical oxygen. Concerning 
the filter layers, the bottom of the BiO-like conduction 
band (or bismuth pocket) along the M(tt, 0) direction 
carries the character of the horizontal p-orbitals of the 
surface oxygens O(Bi) (see Fig. [2] (a)). 

In tunneling calculations, we directly evaluate the 
Green's function instead of diagonalizing the Hamilto- 
nian. For this purpose, the normal state Green's func- 
tion is solved first by starting from the orbital matrix 
elements of the Green's function: 



££(e)' 



(2) 



where e a is the onsite energy of the orbital a. At this 
point, a diagonal self-energy = T, a ' ± iT, a n can be 
included straightforwardly. The simplest self-energy is 
a constant broadening of the states in the form of a 
convergence factor = ^irj. Appendix A (Eq. ( |A2| )) 
presents a more general self-energy which we use to model 
electron-boson coupling. 

The total Green's function G is constructed by solving 
Dyson's equation 

G = g + gVG, 

where V a p are the off-diagonal overlap integrals of Eq. 
([!]) Dyson's equation is exactly solved using the method 
described in Ref. 28 , which is suitable for tunneling 
calculation^. 



C. Pairing interaction and the superconducting 
state Hamiltonian 



Superconductivity is included by adding a pairing in- 
teraction term A in the Hamiltonian of Eq. as follows 



H = Hi 



a(3a 



(3) 



A gap parameter value of |A| = 0.045eV is chosen to 
model a typical experimental spectrum 2 for the illustra- 
tive purposes of this study. We take A to be non-zero 
only between d x 2_ y 2 orbitals of the nearest neighbor Cu 
atoms, and to possess a d-wave form, i.e., A d ( d ± x ^ = 
+ |A| and ^d(d±y) — ~ where d denotes the d x 2_ y 2 
orbital at a chosen site, and d±x/y the d x 2_ y 2 orbital of 
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the neighboring Cu atom in x/y-direction. In momentum 
space, the corresponding A is given by 



[cos k x a — cos k y a] 



(4) 



where a is the in-plane lattice constant. The pairing in- 
teraction of Eq. (|3| allows electrons of opposite spins to 
combine to produce superconducting pairs such that the 
resulting superconducting gap is zero along the nodal di- 
rections k x = zLZfy, and is maximum along the antinodal 
directions. This choice of pairing interaction follows, e.g., 
the one-band formalism given in Ref. [30j 
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FIG. 3: (color online) Main: LDOS (or the diagonal elements 
Polol of the density matrix; see Appendix B for details) of 
d x 2_ y 2 (green) and d z i (blue) orbitals of Cu. [Note d z 2 curve 
is scaled up by a factor of 4.5 to compare the shapes of the 
two LDOSs.] Oscillations at high positive or negative energies 
( above ±0.5 eV) are artifacts due to the use of a sparse mesh 
of k-points in the computation. Inset: Anomalous density 
matrix term p^g discussed in Appendix B, where a and /3 
denote d x 2_ y 2 orbitals of two neighboring Cu atoms. 

For treating the superconducting case, we employ the 
tensor (Nambu-Gorkov) Green's function Q (see Ref. [31} 
with the corresponding Dyson's equation: 



where 



g = g° + gvg°, 



and V 



(5) 





where G e and G^, denote the Green's functions for the 
electrons and holes, respectively. 

The normal state electron Green function G e can be 
used to derive the hole Green function Gh- It can be 
shown by, e.g., the equation of motion method, that 



It is straightforwardly shown then that 

G e = Gg + FA^Gg 
F = G e AG° h 



(6) 



The quasiparticle Green's function projected onto elec- 
tron degrees of freedom is then written in the form 



G e = G° + G e Z BCS G° where E 



BCS 



AG°At. (7) 



We also need the self-energy term for holes. Since 
the transformation from electron to holes follows that of 
the Green's function, we obtain the general form 

E h a (e) = -K*(~e) = -sL(-e) + <(-e). 

In our particular case, we use a self-energy with an odd 
real part and an ev en im aginary part as discussed in Ap- 
pendix A (see Eq. ( |A2| )) Our self-energy is thus invariant 
under electron-hole transformation. 

Figs. [2|e) and (f) show the modifications of the nor- 
mal state band structure from the introduction of the 
pairing interaction. Only the region within ±500meV of 
the Fermi level is shown in panel (e), as the remainder of 
the bands are unchanged from the normal state results of 
panels (a) and (b). The superconducting state dispersion 
in panels (e) and (f) clearly displays a d-wave gap with a 
maximum in the antinodal region near the M point and 
zero gap along the nodal direction near (7r/2,7r/2). Note 
that both bonding and antibonding VHSs possess gaps of 
similar magnitude. Fig. [2|e) also shows the relative elec- 
tron/hole character of the quasiparticles. As expected, 
the quasiparticles are very distinctly either electron- or 
hole-like almost everywhere except within a very nar- 
row energy range at the top and bottom of the SC gap. 
Fig. [3] further shows that mixing of the electron and 
hole features gives rise to coherence peaks in the LDOS 
of Cu-d x 2_ y 2 and to a lesser extent in the LDOS of Cu- 
d z 2 . The effects of electron- hole mixing are however most 
pronouned in the anomalous matrix element of the quasi- 
particle Green's function (inset to Fig. [3] and Fig. [4|. In 
fact, the off-diagonal matrix element between an up- spin 



electron orbital and a down-spin d x 



hole or- 



Gl8 a (-E) 



bital of two neighboring Cu atoms gives the most impor- 
tant term in the anomalous part of the Green's function. 
This term has d-wave symmetry, which manifests itself 
as a change in sign each time we make a rotation of ^ 
around the central Cu site. In addition to the coherence 
peaks, the anomalous density matrix inherits features 
from the VHSs in the regular part of the density matrix, 
which in view of electron-hole symmetry are reflected on 
both sides of the Fermi energy. Additionally, strong hy- 
bridization between up-spin Cu-d x 2_ y 2 electron orbitals 
and down-spin orbitals of O p x holes (and vice versa) 
takes place as shown in Fig. [4] This term is comparable 
in strength to the Cu-d — Cu-d terms and changes sign 
in rotations of ir for reasons explained in the special case 
(3) of the following paragraph. Fig. [4] also shows a small 
onsite contribution from the up-spin electron and down- 
spin hole of the p x -orbital on the oxygen between two 
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neighboring Cu atoms. It is notable that these matrix el- 
ements strictly follow the d-wave symmetry in rotations 
around the central Cu atom. 
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FIG. 4: (color online) Main: Matrix elements of the anoma- 
lous Green's function for onsite p x -orbital of an intermediate 
oxygen atom (green lines), d x i_ y i orbitals of two neighboring 
Cu atoms (red lines), and between Cu-d x 2_ y 2 and &p x orbital 
of a neighboring oxygen (black lines). Inset: The directional 
dependence of the sign of the off-diagonal element F dp . For 
details see special case (3) in the text. 



These transformation properties follow consistently 
from Eq. §6§. Let us, for example, look at the equa- 
tion in the ^-direction: F a/3 = G eyad A d ^ d±x) G° h ^ d±x)(3 , 
where d is a shorthand notation for d x i_ y i of a chosen 
Cu atom, and d ± x stands for the d x i_ y i orbital of the 
neighboring Cu atom in the positive/negative x-direction 
and consider several specific cases as follows. 

(1) For a = d and (3 = d±x, both G e ^ ad and G® h ( d ± x ^p 
are onsite matrix elements, and thus their sign remains 
invariant when changing from one Cu to another. Hence 
the term A d ( d ± x ^ is decisive, and the sign can change 
only in going from x- to y- direction; 

(2) For a = f3 = 0-p x1 we have to first look at the 
term G e , Vxd G\ 

(d±x)p x ' Since the relative phases of the 
off-diagonal matrix elements of the Green's function are 
proportional to the sign of the overlap of the two or- 
bitals, it is straightforward to see from the signs of the 
lobes of the d and p orbitals that this product is invari- 
ant to change in direction as well as in going from x to 
y. Therefore, A d ( d ± x \ again gives the d-wave symmetry 
of these terms; 

(3) For a = d and j3 = 0-p x: G e , ad is diagonal and 
thus invariant. Considering the overlaps, one sees that 



u h,(d- 



— _/^° 

— ^ r h,(d+x)pi 



and 



But, since A d(d±x) = -A d(d±y) , 



Px(-) 



r d Px {+) ~ ~r dpyW - r dVy (_), 

as sho wn in the inset to Fig. [4j 

Eq. ( |B5| ) of Appendix B shows that F a p oc (c a ^cpi). 
Hence, case (3) of the last paragraph indicates that there 
is a significant pairing {c dx2 _ 2^ c <t>i) wnen 

\4>) OC \ Px (+)) + \Py(+)) ~ M")> " K(-)>- 

Recall that we introduced superconductivity in Hamilto- 
nian of Eq. ([3| only on the Cu-d x 2_ y 2 orbitals. Thus 
we see that within our model the strong Cu-0 hybridiza- 
tion automatically induces pairing on the oxygen orbitals. 
This pairing is analogous to the concept of Zhang-Rice 
singlets (ZRS) in the low doping limits 2 -!, where pair states 

t)|<H) - \d x 2_ y 2 t) 

are formed. Note, however, that ZRS is a concept re- 
lated to doping levels in the 'normal' phase, and is not 
directly concerned with superconductivity. Nevertheless, 
the preceding considerations indicate that our model is 
in accord with the ZRS scenario of the normal stated 



D. Green's function formulation of tunneling 
current 

We turn now to consider the formulation of the tun- 
neling spectrum. For this purpose, we apply the c on- 
ventional form of the Todorov-Pendry expressio n^ 1 * 12 * for 
the differential conductance a between orbitals of the tip 
(t, t') and the sample (s, s'), which in our case is straight- 
forwardly shown to yield 



dl 
dV 



2^e 2 

n 



Pw(E F )Vt>sP B AEF+eV)V} v (8) 



tt'ss' 



where the density matrix 
Pss> = -z Im l G ts>] 



2iri 



(9) 



is given in terms of the retarded electron Green function 
or propagator G~^ 3 ,. Eq. ^ differs from the more com- 
monly used Tersoff-Hamann approach^ in that it takes 
into account the details of the symmetry of the tip or- 
bitals and how these orbitals overlap with the surface 
orbitals. 

Since electrons are not eigenparticles in the presence of 
the pairing term, Dyson's equation needs to be applied 
to the Green's function tensor: 

g- = g+ + £ + ($r - = g + - 2ig + ^'g~ (10) 



After extracting the electron part from Eq. (10) and 
applying Eq. ([9|, the spectral function can be written 



U h,(d+y), 



~ U h,(d- 



^h,{d+x)P' 



(11) 



7 



Using Eq. (11), the tunneling current of Eq. ([8| can be where 
recast into the form 



* = E r * 



(12) 



2e 2 



t'ss' 



p u >{E F )V t , s {Gt^\G- s> + F+J?' a F- f )v} t , 



(13) 



r 



and the Green's function and the self-energy are evalu- 
ated at energy E = Ep + eVb. Eqs. (12) and (13) are an 



extension of the Landauer-Biittiker formula for tunneling 
across nanostructures (see, e.g., Ref . l34|) . and represent 
a reformulation of Refs. I35landl36l By comparing Eqs. 



(11) and (13), we see that if the tip makes contact with 
only a single surface atom orbital, e.g., a Bi-p z orbital, 
then the tunneling current is directly proportional to the 
LDOS of that orbital In particular, the tunneling current 
bears in general no such simple relationship to the quan- 
tity of most interest, namely, the LDOS on the Cu02 
plane. Obviously, the tunneling formalism of Eq. (13) 



must be further elaborated in order to find the relation 
between the interesting LDOSs and the tunneling spec- 
trum. 



1. Tunneling channels, filter function and tunneling matrix 
element 

The experimental STM spectra in the cuprates have 
to date been mostly compared to the electronic LDOS of 
the superconducting cuprate layer, especially the LDOS 
of the Cu-d x 2_ y 2 orbital. The discrepancies between the 
spectra and the LDOS are then ascribed to 'tunneling 



J 



matrix elements' or 'filtering functions :8 . The former 
refers to the general problem of modeling spectroscopies, 
where the signal is distorted by the spectroscopic pro- 
cess, and may even vanish due to the presence of selec- 
tion rules. The latter term refers to how the states of 
electrons (or quasiparticles) from the initial state within 
the superconducting layers are modified when traveling 
through the oxide overlayers before reaching the tip. Eq. 
( 13 ) above accounts fully for the tunneling process, and 



it can be reformulated to reveal, for example, the filter- 
ing effect more clearly. For this purpose, it is convenient 
to the denote various orbit als as follows: s and s f for the 
orbit als of the sample surface, which overlap with the tip 
orbital t\ f and /' for the orbitals of the filter layers, BiO 
and SrO; c and c' for orbitals in the cuprate layer; and, 
a for any orbital that is singled out, which in our case 
usually will be an orbital in the cuprate layer. Denoting 
the Green's function for the filter layers decoupled from 
the rest of the system by G^J , and the matrix elements 
within the cuprate layer in the coupled system by G+,, 
application of Dyson's equation to Gf a yields 



Gta = G°+V fc Gt a and F+ = G°+V fc F ( 



+ 

COL 



Hence, Eq. (13) can be written as 



2e 2 



Ptt>(E F )M t , c (Gt a X"aG- c , + F+j:" a F- c ,)Ml t 



(14) 



t'cc' 



r 



where 



M tc = V ts G°+V fc , 



(15) 



which gives the filtering amplitude between the cuprate 
layer and the tip, and constitutes a multiband generaliza- 
tion of filtering function of Ref. El Similarly, the matrix 
element of the density of states operator p cc > within the 
cuprate plane can be recovered in terms of the spectral 



function: 

27re 2 



h 



Ptt>(E F )M t , cPcc ,(E F + eV)M\, v (16) 



Eqs. (14)-(16) show a number of interesting aspects of 



the tunneling process as follows. 

(1) Since applying the filtering matrix element M tci 
which describes the effect of the BiO and SrO overlay- 
ers, involves M and M\ interference effects will occur 
between various paths to the tip from the cuprate layers 
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through the filter layer; 

(2) The partial current terms in ( 16 ) under the summa- 
tion are proportional to elements of the density matrix 
confined to the cuprate layer. Only orbit als with a no- 
table overlap with the p z orbital of the apical oxygen on 
the SrO layer will give a significant contribution to the 
total current; 

(3) The partial elements of the spectral function 



(a) Experiment (b) Corrugation 



:( G tc^"aG ac , + F+ a Z" a F ac ,) 



(17) 



extracted from Eq. ( 14) show which orbitals a contribute 



to the chosen element of the density matrix p cc > . Further- 
more, the current contribution T ta between the tip can 
be divided into regular and anomalous terms and , 
respective!}^. 

Since the filter layers are insulating at low energies, 
these layers will give little structure to the spectrum at 
low bias voltages, so that the structure of the spectrum 
is mainly controlled by the matrix elements p cc /, and in 
this sense the spectrum is a filtered mapping of the LDOS 
of the cuprate orbitals. We will show however that the 
Cu-d x 2_ y 2 orbitals right below the tip do not enter the 
spectrum through Eq. (16) since their overlap with the 



relevant orbitals of the SrO layer is zero. Instead, Cu-d z 2 
has a large overlap with p z of the apical oxygen and hence 
these orbitals of the Cu atoms play a dominant role in 
the tunneling spectrum. 

The detailed contribution of any specific orbital a can 
be extracted from Eq. (14). The regular and anomalous 



matrix elements of the spectral function, G+ a Yi" a G ac , 
and F c + Y," a F~ c ,, describe propagation of electrons or 
holes within the cuprate layer from orbital a to the or- 
bitals c and c' . The latter orbitals act as "gates" between 
the cuprate layer and the filter layer. For example, if a is 
d x 2_ y 2 of a Cu atom and c and c' are d z 2 orbitals, which 
strongly overlap with the filter layer, the matrix element 
filtered by M and gives the contribution of a specific 
d x 2_ y 2 orbital to the total tunneling spectrum. Note that 
in the superconducting state the anomalous matrix ele- 
ments of the spectral function must also be considered. 
Fafiij) involves the creation of an electron with spin up 
coupled to the annihilation of a hole with spin down given 
by (ct | (r)c^(O)), and thus describes the formation and 
breakup of Cooper pairs as shown in Appendix B. The 
decomposition of Eqs. (14)-(16) are, in fact, a general- 
ization of the tunneling channel approach to transport 
through one- molecule elect ronic components^ and STM 
of adsorbate molecule J 39 | 4Q | . In the present context, the 
"tunneling path" analysis gives us the "origin" of the 
signal, since G^T/' a G~ c , gives the probability of propa- 
gation between orbitals a and c. 



60 A 




FIG. 5: (a) Typical experimental topographical STM map 
after Ref. [4] (b) The computed corrugation of two STM line 
scans and (c) theoretically predicted topographic map. The 
two paths are shown in (c) by arrows. 



III. RESULTS 

A. Topographic maps 

We discuss first the topographic STM map, i.e., the 
constant current surface for a tip scanning across the 
sample surface. The computed topographic map is very 
robust against changes in measuring parameters such as 
the bias voltage or the tip-surface distance. Figure [5] 
compares the calculated and typical experimental results. 
Furthermore, corrugation along two paths of line scan is 
shown in Fig. [5] (b). The Bi atoms are seen as bright 
spots, while the surface oxygens are dark due to very low 
current coming through these surface atoms. We will see 
in connection with the analysis of the tunneling channels 
below that the apical oxygens act as the primary gate for 
passing electrons from the Cu02 layers up to the surface 
BiO layer. Accordingly, the Bi atoms appear bright be- 
cause there exists an easy channel between the surface 
Bi atoms and the apical oxygens below via the Bi p z or- 
bitals. On the other hand, the oxygens in the surface 
layer are dark because the p XiV orbitals of O(Bi) are or- 
thogonal to the (assumed) s-symmetry of the tip, while 
the O(Bi) p z orbitals are relatively weakly coupled to the 
p z of the apical oxygen as discussed below in connection 
with Fig. [f| 



B. Tunneling spectra 

Fig. [6] (a) compares a typical experimental (red line) 
STS spectrunP to the calculated one (black line). The 
overall agreement between theory and experiment is seen 
to be good, although the VHSs are seen as separate struc- 
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tures in the calculated curv d 41 * 4 ^ . The agreement also 
extends to the low energy region shown in Fig. (6^b), 
where the width and positions of the coherence peaks is 
reproduced reasonably well. 43 The tendancy for increas- 
ing intensity towards negative bias is seen in both mea- 
surements and computations. This is in sharp contrast 
to the shape of the LDOS of Cu-d x 2_ y 2 orbital (green 
curve). As emphasized in Ref. 10 , this remarkable asym- 
metry of the spectrum between positive and negative bias 
voltages reflects the opening up of channels other than 
Cu-d x 2_ y 2 , especially of Cu-d z 2 , as one goes to high nega- 
tive bias. This asymmetry thus appears naturally within 
our conventional picture and cannot be taken to be a hall- 
mark of strong correlation effects as has been thought to 
be the case. 
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FIG. 6: (color online) (a) A typical experimental tunneling 
spectrum (red line) from Bi2212 (after Ref. [2} is compared 
with the calculated spectrum (black) . The green curve shows 
the LDOS of the Cu-d x 2_ y 2. (b) Expanded view of the ex- 
perimental and calculated spectrum in the low ener gy region, 
(c) Comparison of the model self-energy (Eq. (A2)) assumed 
for the C\x-d x 2_ y 2 orbitals and the self-energy from the con- 
volution of a De bye-t ype phonon spectrum and the LDOS of 
C\x-d x 2_ y 2 (Eq. (Al)) as discussed in Appendix A. 



There has been considerable interest in understanding 



the coupling of electrons to bosonic modes in the so-called 
'low-energy kink' region within ~ ±100 meV of the Fermi 
level. In particular, the peak-dip-hump structure seen in 
the experimental spectrum in Fig. [6^b) is generally be- 
lieved to be the result of the coupling of electronic degrees 
of freedom to a collective mode (Refs. I9)44|45p . Fig. 
|6jb) shows that the peak-dip-hump feature can be de- 
scribed by our simple self-energy correction discussed in 
Appendix A. This point however requires further study, 
including an analysis of how this feature evolves with 
doping. 



C. Selection rules 

The filter function M tc controls selection rules dictated 
by matching of the symmetry properties of the cuprate 
layer, filter layers and the tip. A closer examination of 
M tc reveals that strong tunneling through the apical oxy- 
gen layer is associated with a matching of the symmetry 
of the cuprate layer wave function to that of the api- 
cal 0-p z . The key is the relative symmetry of the wave 
functions with respect to the axis of tunneling: An 'odd' 
wave function, e.g., the Cu-d x 2_ y 2 has zero overlap with 
an 'even' wave function such as 0-p z . In contrast, two 
orbitals with the same symmetry couple more strongly. 
Accordingly, the p z of the apical oxygen and the Cu-d z 2 
possess large overlap, while Cu-d x 2_ y 2 has zero overlap 
with any s- or p-orbital of the apical oxygen. This is the 
reason that direct tunneling is forbidden between Cu- 
d x 2_ y 2 and the s-wave symmetric tip through the filter 
layer. Hence, M tc functions here are consistent with the 
filter function of Ref. [HJ Similarly, coupling between an 
s-wave tip and the p x and p y orbitals of the Bi atom lying 
directly below the tip is forbidden. Therefore, within the 
filter layer, the main 'vertical' overlap is between the p z 
orbitals of Bi and apical oxygen, and these orbitals indeed 
are found to provide the main channel through the filter 
layers as depicted in Fig. [7fa). We find additional rel- 
atively small contributions from the on-site Bi-s-orbital 
and p-orbitals of the surrounding Bi and O(Bi) atoms, 
but such 'background' contributions to the current do 
not seem to be dominated by any particular channel. 

Figure [7Jb) illustrates another example of a symmetry- 
forbidden tunneling path, where the tip is centered be- 
tween two surface Bi's, i.e. on the top of an oxygen of 
the cuprate layer. Since we assume an s-wave tip with 
negative hopping integrals to the nearby Bi atoms, when 
we follow either path up to the Cu-d z 2 orbitals, the signs 
of the hopping integrals are identical. However, the 0-p x 
orbital between the two Cu atoms changes sign from one 
Cu to the other. This gives the two paths from 0-p x to 
the s-wave tip an opposite phase leading to destructive 
interference between the paths, making the O atom invis- 
ible. However, if the s-wave tip is replaced by one with, 
e.g, p x symmetry, the oxygen would become visible and a 
weaker signal would appear from the neighboring d x 2_ y 2 
orbitals. Experimentally, this could be accomplished by 
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functionalizing the tip by attaching a suitable molecule 
to the tip. A similar procedure has been used to obtain 
a cont rast i nversion for CO molecules adsorbed on a Cu 
surfactPHa 




(b) 



FIG. 7: (color online) (a) Dominant tunneling channel from 
the cuprate layer, from Cu d x 2_ y 2 orbitals through the neigh- 
boring Cu d z 2 to Bi p x to the tip. (b) An oxygen atom in the 
cuprate layer is invisible to a STM tip right above, since the 
paths through Cui and Cu2 interfere destructively. 



D. Tunneling channels 

The origin of the current from the cuprate layer can 
be understood by inspecting the individual terms of Eq. 
(17), which we refer to as 'tunneling channels', i.e., from 

and 



the regular and anomalous elements G^Y*" a G Q 
F^Y," a F~ c , , of the Green's function. [Although tunnel- 
ing channels are a normal state property, the anomalous 
matrix elements play an important role in generating the 
coherence peaks and thus are relevant more generally.] 
For simplicity, we assume that the tip is right above a Bi 
atom. The dominant element of the filter function M t c 
is then between the tip orbital and the d z i orbital of the 
upper layer Cu atom lying beneath the surface Bi atom, 
so we take c — d = Cu-d z 2 in results shown in Figs. [8] 
and [9j Fig. [8] shows the relative contributions of the 
regular and anomalous matrix elements. The near Fermi 
energy current is primarily associated with the d x i_ y i 
matrix elements. While the regular matrix elements of 
Cu-d x 2_ y 2 are almost solely responsible for the spectrum 
at energies around the VHSs, the anomalous elements de- 
termine the features around the gap region, especially the 
coherence peaks. Fig. [8] shows that coherence peaks are 
inherited from the anomalous and not the regular part 
of the Green's function, and reflect physically the effects 
of non- conservation of the number of electrons near the 
gap region. 
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FIG. 8: (color online) Partial spectrum with c — c — d z i in 
Mtc- The regular (red line) and anomalous (green line) com- 
ponents are shown together with the total contribution of the 
two parts (solid black). Blue curve shows the corresponding 
regular Cu-d x 2_ 2 contributions. 



In Fig. [9j the current of d x 2_ y 2 character is further 
broken down into contributions from various neighbors of 
the central Cu atom of the first and second Cu02 layer 
away from the free surface. We see in panel (a) that the 
upper Cu02 layer is more important than the lower one, 
but that the upper layer is by no means dominant. It 
seems that the coupling between the tip and the lower 
layer is strengthened via the relatively large overlap be- 
tween the d z 2 orbitals of the central Cu atoms of the two 
layers, which opens an important interlayer channel. The 
d x 2_ y 2 orbitals of the two layers mix not only to induce 
the well-known bilayer splitting in Bi2212, but also play 
a significant role in the flow of current to the tip from 
the lower cuprate layer. 

It can be seen from Fig. [9] (a) that the d x 2_ y 2 orbitals 
of the four nearest-neighbor Cu atoms of the central Cu 
give a significant contribution to the total spectrum, but 
that this amounts to only about one third of the con- 
tribution from all d x 2_ y 2 terms from the upper layer. 
Due to the non-local nature of the Bloch-states within 
the cuprate layers, it is clear then that the total sig- 
nal involves long range contributions, and attributing the 
spectrum merely to the four nearest neighbor Cu atoms 
provides only a rough approximation. 

Anomalous contributions are considered in Fig. [9] (b) . 
Here, the upper and lower layers give an almost equally 
large contribution, indicating that coherence peaks also 
are not all that local in character. Notably, we find a fi- 
nite onsite anomalous contribution of d x 2_ y 2 even though 
the regular term is zero. This can be understood with 
reference to Eq. Consider the term 

F z 2 d = G° Z 2^ d+Xi )^(d+ Xi )dGh,dd^ 



where d is shorthand for d T 



of the central Cu and 
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FIG. 9: (color online) (a) Various contributions to tunneling 
spectrum from the regular matrix elements (assuming c — 
c' — d z 2 in Mtc), Tfo, of Cu-d x 2_ y 2 orbitals of upper and 
lower Cu02 layer. Contributions from the nearest neighbor 
(nn) Cu atoms in the upper and lower layer are shown, (b) 
Same as (a), except this panel refers to the contributions from 
the anomalous matrix elements, T^. 



representations of the two-dimensional C± v group, the d- 
orbitals \d x 2_ y 2,i) of the site i participate in eigenfunc- 
tions of the system as a linear combination 



\d x 2_ y 2,i). 



This combination of the four neighboring orbitals at 
(0, ±7r) and (±7r, 0) belongs to the same representation 
of C± v as the 4s and d z i orbitals of the central Cu atom 
(see Fig. [To]), as well as the p z orbitals of the apical 
oxygen and the surface Bi atom. At this k-point, the 
phase difference between the lattice sites causes all the 
d-orbital lobes pointing towards the central atom to have 
the same sign. Hence, this combination yields a large off- 
diagonal element overlap with the surface p^-orbital, and 
a dominant tunneling contribution around the gap. Simi- 
lar arguments can be applied to understand contributions 
from other farther out atoms. An example was given in 
Fig. [7^b) above where the position of the tip and the 
symmetry of the relevant orbital strongly influence the 
visibility of an atom. 



d + Xi is d x 2_ y 2 of the neighboring Cu in either x- or y- 
direction. Clearly, G Q e z 2^ d+x .^ transforms under rotations 
of | in the same way as A( d+X .) d , and since G^dd 1S an 
onsite term, the combination is invariant. Hence the four 
terms in the sum over the neighbors are equal, yielding 
a non-zero onsite term. 

We emphasize that the anomalous contribution of the 
four neighboring Cu atoms is quite small. Let us consider 
the term 



F. 



z 2 {d+ Xi ) 



— ^e,z 2 



l^d(d+Xi)Gh,(d+Xi)(d+Xi)- 



Due to symmetry, G° e z2d = 0, and thus this term van- 
ishes. However, there are terms like 

Fz 2 (d+ Xi ) = Ge i z 2 (d+2xi)A(d+2x i )(d+Xi)Gh,(d+Xi)(d+Xi) 

which do not vanish, but are very small, since G° e z 2^+2x-) 
is a relatively small term. A similar analysis can be car- 
ried out for the second and third neighbors. The second 
nearest neighbors, which lie along the nodal direction 
in k-space, give the largest single contribution, although 
this contribution is not dominant. The third neighbor 
contribution is a little larger than the onsite contribu- 
tion. 



IV. FURTHER COMMENTS 

A. Symmetry Analysis 

The selection rules can be formalized using group the- 
oretical arguments related to the filtering functionP. For 
example, in order to explain the dominance of the d x 2_ y 2 
orbitals of the four neighboring Cu atoms, considering 



(a) T-point (0,0) (b) M-point (0,7i/a) 

^jLg ^? 2 opposite phase ^^j[jf" y2 
W s ame phas e w 

+ + 



FIG. 10: (color online) Relative phases of the central d z 2 
orbital and the neighboring d x 2_ y 2 orbitals at the F point (a) 
and at the M point (b). 



B. Electron extraction/injection 

To relate the tunneling current to the LDOS of the 
cuprate layer, we have introduced the concept of tun- 
neling paths through Eq. (14), which implies that each 



path begins or ends on a particular atomic orbital. This 
non-intuitive concept requires some comment. In reality, 
the current flows through the sample with each electron 
ejected to the tip being replaced by an electron from a 
distant counterelectrode. For a simple system, such as 
a nanostructure, non-equilibrium Green's function for- 
malism with two 'leads' closing a current circuit have 
been invoked (see, e.g., Ref. l34|) . Tersoff-Hamann (TH) 
or Todorov-Pendry (TP) approach, on the other hand, 
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assumes that the current is composed of a series of tun- 
neling events^, and that the replacement of electrons at 
the counterelectrode has a negligible effect on the tun- 
neling process. Since the current in STS is of the or- 
der of 10 — lOOpA, there is only about one electron each 
1 — 10ns which flows across the sample, justifying the as- 
sumptions underlying TH/TP approach. Both TH and 
TP are based on calculating individual tunneling events 
in a LEED-like formalism 49 . Due to the finite an 
electron created on a particular atom will have only a fi- 
nite probability of escaping to the tunneling tip, and Eq. 
(14) shows how to add up the contribution of all these 



tunneling processes in terms of the equilibrium LDOS of 
the sample. 



V. CONCLUSIONS 

We have presented a comprehensive framework for 
modeling the STS spectra from the normal as well as the 
superconducting state of complex materials in a material- 
specific manner. Our formulation makes transparent the 
connection between the LDOS and the STS spectrum or 
the nature of the tunneling 'matrix element', and it is 
cast in a form that reveals the filtering effect of the over- 
layers separating the tip and the layers of interest. Our 
decomposition of the tunneling current into contributions 
from individual local orbitals allows us to identify impor- 
tant 'tunneling channels' or paths through which current 
reaches the STM tip in the system. Our analysis high- 
lights the importance of anomalous terms of the Green's 
function, which account for the formation and breaking 
up of Cooper pairs, and how such terms affect the STS 
spectrum. 

We apply the formalism to the specific case of Bi2212. 
Mismatch of symmetry between orbitals on adjacent 
atoms, or between the tip and the sample orbitals, is 
shown to severely restrict the corresponding contribu- 
tion to the tunneling current. For these reasons, the 
contribution from Cu-d x 2_ y 2 orbitals comes not directly 
from the Cu-atom lying right below the Bi atom, but 
from a fourfold symmetric indirect route involving the 
four nearest-neighbors of the central Cu as well as longer 
range background from farther out Cu-d x 2_ y 2 orbitals. 
In the superconducting state, the coherence peaks of the 
spectrum are shown to be dominated by the anomalous 
spectral terms, which also are found not to be all that 
localized around the central Cu atom. In particular, we 
find a small anomalous on-site term and a practically 
vanishing first nearest neighbor contribution, with most 
of the anomalous contribution arising from the second 
neighbors and beyond. 

We have concentrated in this study on the large hole 
doping regime of the cuprates where a homogeneous elec- 
tronic Fermi liquid phase is consistent with most exper- 
iments. The fact that we have obtained good overall 
agreement between our computations and the measure- 
ments, especially with respect to the pronounced asym- 



metry of the spectrum between positive and negative bias 
voltages, indicates that this remarkable asymmetry can 
be understood more or less within our conventional pic- 
ture without the need for invoking exotic mechanisms. 
At lower dopings, strong correlation effects including the 
possible presence of competing orders or inhomogeneous 
electronic states (nanoscale phase separation) would need 
to be taken into account. However, the present frame- 
work can be extended fairly straightforwardly through 
the addition of Hubbard terms in the Hamiltonian to 
provide a viable scheme for investigating the tunneling 
response throughout the phase diagram of the cuprates 
and other complex materials, including the modeling of 
effects of impurities and dopant atoms in the system. 
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Appendix A: Boson-electron coupling 

In the vicinity of the Fermi energy, dispersion anoma- 
lies are found in ARPES spectra arising from coupling 
of electronic degrees of freedom to phonons and/or mag- 
netic modes, often giving the appearance of a peak-dip- 
hump feature^. These boson-electron couplings also 
strongly affect the STS spectrum^. This appendix dis- 
cusses a model self-energy for describing such anomalies. 

A significant contribution to the electron-phonon cou- 
pling is associated with modulation of the electronic hop- 
ping integrals by the phonons. The generalized coordi- 
nate of atomic displacement in g-basis is quantized in the 
standard way: 



Qq 




where a q (a q ) is the annihilation (creation) operator of the 
phonon mode q, and Q q is the frequency of the mode. 
However, the most natural way to couple this to real- 
space tight-binding basis is to make a transformation to 
the basis of real space displacement of atom a in the 
following way: 

<V = (v\q)Q q , 

where Einstein summation over phonon modes q is im- 
plicit. Note, that \i is a composite index denoting both 
the index of an atom and the direction of displacement. 
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Consequently, in tight-binding basis, this gives rise to 
a term in the Hamiltonian of the form 



H' 



el—vib 



dV a6 A t 



where V a s is the hopping integral between orbitals a and 
is the coordinate of atom ja. 
This coupling can be embedded into the electronic 
Hamiltonian as an energy dependent self-energy. Fol- 
lowing the arguments of Ref. [50j the general form of 
self-energy is written as: 



((1 - f(e - m) + n b (n))G±(e - hfl) 

+(/(£ + nn) + n b (a))G%(e + nn)), 



(Al) 



where g^ v (P) = ^ q (^\q)5(Q — Q q )(q\v) is an element of 
the vibration mode density matrix. Note again that we 
use Einstein summation convention, so that summation 
is implied over orbital indices 7 and 5 and the phonon 
polarization indices p and v. 

For simplicity, we now assume that: (i) The bosonic 
coupling only affects the Gw-d x i_ y i orbitals, where we 
include a diagonal self-energy of the form, g(fl) = g£l 2 
when Q < and it is when Q > Q^. For a De- 
bye spectrum of phonons, Qd is the Debye cut-off fre- 
quency, and the normalization factor is g = (ii) 
ps 1 = — ^Im[G^] = pis approximately a constant. This 
amounts to assuming that the electronic density of states 
is smoothly varying within the range of the phononic 
spectrum; (iii) Take T^ 5 = T, a constant parameter. Us- 
ing these assumptions, the final form for the self-energy 
is 



E+ 



■£(<a. + *r> + («"-i)htei 



(A2) 



where z = (e + irj)/(fiQd), A = Jf^T 2 /), and 77 is a 
convergence parameter. Although we have derived the 
preceding form for coupling to a 3D Debye spectrum of 
phonons, the results are not too sensitive to details of 
the spectrum, and we would expect a similar result for 
an Einstein phonon or the magnetic resonance mode 5 -^. 

It is interesting to consider the asymptotic forms of 
self-energy as follows. If fiVtd <^ 



A I — 

7TZ 



For large boson energies, i.e., ftVtd ^> £> we obtain 
E(e) « -A l-z + iz* 

\ 7T 



(A3) 



While Eq. (Al) gives a general form of phononic self- 



energy for any pair of orbitals, in the present calculations, 
we adopt a few simplifications. First, we assume only 



diagonal terms of self-energy to make the model more 
tractable. Second, we apply Eq. (A2) to Cu-d x 2_ y 2 or- 
bitals using parameters TiVtd — SOmeV and A = 60meV. 
The former value gives the best fit to the peak-dip-hump 
structure, and the latter controls the smoothness of the 
spectrum. In Fig. j6^c) we make a comparison between 
the more general form of Eq. ( Al ) with the accurate den- 



sity of states of Cu-d x 2_ y 2 orbitals. For the remaining 
orbitals we mimic a Fermi-liquid type self-energy, which 
can be modeled with a E" ex e 2 and E 7 oc e; here we 



employ the asymptotic form of Eq. (A3), choosing pa- 
rameters hQd = 2.0eV (to ensure the correct asymptotic 
form for whole the energy range) and A = lOOmeV. In 
this way, the need for a Kramers-Kronig transformation 
is avoided. 

We can straightforwardly include in the self-energy the 
effect of magnon scattering 42 responsible for the high en- 
ergy kink 52 . This will broaden the spectrum in the vicin- 
ity of the VHS peaks, thereby improving agreement with 
experiment in Fig. [6^a). It should be noted, however, 
that a more accurate modeling of the self-energy will be 
required both for the bosonic coupling and the Fermi- 
liquid term for treating the underdoped system. 



Appendix B: Bogoliubov quasiparticles in 
tight-binding basis 

This appendix discusses aspects of the Bogoliubov 
transformation within a tight-binding basis. The Bogoli- 
ubov transformation is not explicitly carried out in the 
present calculations since the Green's function tensor is 
obtained directly from Dyson's equation. Nevertheless, 
understanding the relation between the transformation 
and the Green's function tensor in the tight-binding ba- 
sis is necessary for interpreting some of our results. In 
particular, our analysis of pairing symmetry is based on 
the relation between F a p and (c a ^cm). 

The Bogoliubov transformation 53 is conventionally 
carried out in a combined basis of spin-up electrons and 
spin-down holes: 



Cfet 



(Bl) 



-kl' 



These c's diagonalize the one-particle Hamiltonian of Eq. 
([!]) via the transformations 



and 



C U = <-%> c -n = ( a \ k ) c -kV 
or in a more compact form: 



Cat 



(a\ 



(a\ 



Cfc = i^fcC/c, (B2) 
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with inverse c k = B ka c a . 

This change of basis diagonalizes the one-particle 
Hamiltonian: 

e k = (k\a)H liaP (fi\k) 

(with summation over a and f3). In this basis the Hamil- 
tonian of Eq. ([3| becomes 

H = e fc 4 t c fet -h^(l-c_ H cl H )+Afc4t c -H+ A l c -H c fct 5 

now with summation over k. After shifting this by a con- 
stant energy, it assumes the simple form 



H eff = c t jj 



where 



H 



.Al -e k )' 
which can be diagonalized into 

H eff = c *U- 1 UHU- 1 Uc, 

where 



(B3) 



U 



U * Vk and U- 
-v* k u k I 



1 = / Uk ~Vk 



The coefficients are chosen in the standard way in order 
to obtain a diagonal matrix 

UHU- 1 = ( Ek ) , 

o -E k r 



with E k = y^l+jAkf. 

This Bogoliubov transformation introduces the quasi- 
particle basis 



Uc. 



Since we are working in the tight-binding basis, we end 
up with 

( a k \ = ( u* k (k\a) v k (k\/3) \ f c at \ 
(summation over a and (3) or inversely 



V4y {(mm (m< ) \ bi k ) 

(summation over k). 

We are particularly interested in writing the expecta- 
tion values of electron and hole densities, (c^c^), and 
(caacL), and pairing amplitudes (c^cj^), and (c^c at ) 



in terms of the Green's function tensor. For this purpose, 
we start with a 2 x 2 tensor 

(c a cj) = (B ak c k clB k p) = (BakU^asJUBw). (B4) 

Using the fact that (a k a k ) = 1 — f(E k ) and (b\b k ) = 
f(E k ), we evaluate each element of the tensor (c a c^) 
separately as follows: 
(1) The number density 

<4t c /3t) = <0I*> (M 2 /(£fc) + M 2 (l - f(E k ))) (k\a) 
Now we use a trick following Ref. [54] where 

(f3\k)\u\ 2 f(E k )(k\a) = J def(e)(f3\k)u5(s-E k )u*(k\a) 

and 

5{e - E k ) « --Im( ). 



7r e - E k +ir] 



Hence 



where 



(^|fc)|«|V(£?fc)(*l«> = / def(e)p 0a (e), 



p% a {e) = --Im(G+ pa (e)), 

where refers to the electron part of the Green's 

function, 



(a,e\k)(k\e,P) (a\k)\u k \ 2 (k\/3) 



1/3 s-E k + irj e-E k + irj 

It is straightforward to show that 

<4|W = / de\f{e)p% a {e) + (1 - f(e))p h pa (e)], 

where p^p is the hole density matrix. The first part of 
the integral, in fact, gives the number of electrons with a 
chosen spin. The latter part gives the same result as the 
former since the Bogoliubov transformation reflects the 
electron bands to hole bands with respect to the Fermi 
energy; 

(2) The pairing amplitude 

(com) = (a\k) (u(f(E k ) - (1 - f{E k ))v) (k\/3). 

Using the trick of Ref. 54 again gives us the formula 

(com) = - J de{\ - 2f(e))p%(e), (B5) 



where 
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and 



KM 



(a,e\k)(k\h,(3) _ (a\k)u k v k (k\p) 



af3 ^ / s-E k + ir] e-E k + ir] 

In the same manner, one can see that 

Mi) = - / <fe(l - 2/(e))pt f (e), (B6) 



where 



pt^) = - l -Im{{F% a {e)), 



Equations (B5) and (B6) also reveal how the anoma- 



lous part of the Green's function tensor is related to the 
pairing amplitude (c a ^cpi) in a tight-binding basis, or 
equivalently how the anomalous part of the current is re- 
lated to the making and breaking of Cooper pairs. In 
particular, symmetry properties of F a p are seen to be 
related directly to those of (c a ^cpi). 
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